Phonon momentum and damping of mechanical resonators 
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The concept of physical momentum associated to phonons in a crystal, complemented with some 
fundamental reasoning, implies measurable effects in crystals even at a macroscopic scale. We 
show that, in close analogy with the transfer of momentum in the kinetic theory of gases, physical 
momentum carried by of phonons couples the thermal and the velocity field in a vibrating crystal. 
Therefore an heat flow applied to a vibrating crystal can sustain or damp the oscillation, depending 
on the interplay between the temperature and the velocity gradient. We derive the general equations 
of this effect and show that its experimental confirmation is within reach of current technology. 
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In a crystalline solid, a phonon of wave vector k carries 
a " crystal momentum" ?ik that participates in the inter- 
action with particles or other elementary excitations that 
similarly carry their own crystal momentum. In fact the 
ions of the crystal can be seen as scattering centers, and 
the waves emitted from each ion have constructive inter- 
ference only if the total crystal momentum is conserved 
up to an arbitrary reciprocal lattice vector. Therefore 
crystal momentum is not related to a physical momen- 
tum eventually associated to phonons, as they describe 
the relative internal motion of the atoms, and are not 
associated with any linear momentum of the crystal as 
a whole. The limits of this description becomes evident 
when the translational invariance of the crystal is bro- 
ken, as in a finite-size body, or in case of randomly dis- 
tributed scattering centers. In this case the anharmonic 
terms couple the relative internal motion with the center 
of mass motion, and phonons carry both crystal momen- 
tum Sk and the physical linear momentum p = 37?ik 
[ll-Q. The connection between these two forms of mo- 
mentum is a measure of the anharmonicity and is given 
by the Griineisen parameter 7 of the lattice, which mea- 
sures the relation between the volume of a crystal lattice 
and its vibrational properties. 

The intimate connection between phonon momentum 
and thermal expansion has been first pointed out by Bril- 
louin who showed that thermal expansion can be 
viewed as the reaction of the crystal lattice to the in- 
ternal pressure due to phonons. More recently, an ex- 
periment has demonstrated that the collision of phonons 
can move a small cargo along a carbon nanotube jSj] . The 
concept of corpuscular phonon allows also a theoretical 
description of some dynamical phenomena in thermoe- 
lasticity, such as dynamical thermal expansion in thin 
optical windows Q, and is at the basis of the phonon 
pumping mechanism recently proposed for the cooling 
of nanodevices Here we show that phonon momen- 
tum transfer couples the thermal and the velocity field 
in a vibrating crystal and we predict measurable effects 



in crystals even at a macroscopic scale. 

This issue is relevant for a thorough understanding of 
oscillating systems of every size and for setting the limits 
to the performances of many instruments and devices. In 
fact resonant mechanical systems serve as a basis for a 
number of sensors and fundamental studies (S-llj. Their 
size covers orders of magnitude, from micrometers to 
some meters, and their motion is detected, by electri- 
cal or optical readouts, with high sensitivity, sometimes 
limited only by the uncertainty principle. In these in- 
struments a low damping of the resonant motion is often 
crucial to achieve the desired performances, as in the case 
of gravitational wave detectors [l^ or micro-mechanical 
resonators 13]. Here the laser beam used in optical read- 
out usually heats some parts of the resonator: thermal 
gradients up to 10 K are expected in some components 
of the future cryogenic gravitational wave detectors [lij , 
while gradients of about 1 K are currently applied in 
cryogenic micro- resonators with optical readout Uj . 

As discussed above, a physical linear momentum p = 
37ksfi-A:ik is associated to a phonon of the branch s, 
where ik determines the direction of propagation of the 
phonon and 7ks is the Griineisen parameter of the normal 
mode. As long as we deal with crystals at temperatures 
much smaller than the Debye temperature (for instance 
Qd — 645 K for silicon), we can use the simple dispersion 
relation a;s(k) — ck and write p = 37ksfi'^ik- Accord- 
ingly, an equivalent mass rrik — Sjksfi^, arising from the 
coupling of the oscillatory motion of the atoms with the 
center of mass of the crystal, is associated to a phonon 
moving with velocity cik. In case of a phonon in a mov- 
ing crystal, this relation must be complemented with a 
term depending by the motion of the center of mass. If 
we push forward the analogy between a phonon and a 
particle, it is natural to assign a momentum: 
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P = ^Iksfi^ (cifc + vq) 



(1) 



if the phonon is in a crystal moving with velocity vq in 
respect to an observer at rest. This relation transforms 



the physical momentum between reference systems in rel- 
ative motion, considering that the velocity of the phonon 
relative to the crystal is constant. If the phonon is travel- 
ing in a crystal having parts moving with different veloc- 
ities, we expect that it will couple in sequence with the 
parts of the crystal crossed by its trajectory, changing its 
frequency and/or direction as it moves. As a result of this 
interaction, the phonon will actually transfer momentum 
across the crystal. Even if this issue deserve dedicated 
molecular dynamics studies, as a first approach we as- 
sume that the phonon transfers to the crystal its energy 
and momentum when suffering a collision, starting a new 
path with the frequency and momentum corresponding 
to the local temperature and velocity. 

To understand the physical meaning of these assump- 
tions we propose a simplified model in Figure [1] An 
oscillating cantilever beam is heated from its free end, so 
that both the temperature and the amplitude of the ve- 
locity increase along the x— axis. We consider three thin 
sections of thickness Ax '^If, say Qj with j = 1 — 3, at 
intervals of length If along the cantilever, with If the 
Mean Free Path (MFP) of the phonons. Due to the 
temperature dependence of the phonon density, a larger 
number of phonons and with higher energy are emitted 
from the section in respect to those emitted from Qi. 
Following the dominant phonon approximation |15i] . we 
assume that the whole thermal energy in a section can 
be ascribed to phonons of a single frequency (D. Then, in 
the time dt, we have (ni, ^2, ^3) phonons with frequency 
(wi,(D2,W3) starting respectively from (Qi, Q2, Q3), with 
fis > n2 > fii and oj^ > 0J2 > u)i (Figure [1];) . Now we 
consider the balance of the momentum pz directed along 
the z— axis, at the central section Q2 when the cantilever 
goes through a zero crossing. We assume that all phonons 
starting from Qi and Qa, travel for a length If and then 
suffer a collision with the lattice in Q2. According to Eq. 



([T]), phonons carry a momentum proportional to 



then: 
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where we take in account that only one half of the 
phonons starting in the Qi^Qa are emitted towards Q2. 
It is evident that, in a linear approximation of the ve- 
locity V2 = {vi -\- W3)/2, we have p. = in the isother- 
mal case, when phonon density and frequencies do not 
changes along the cantilever. In this case the excess mo- 
mentum coming from Q3 is compensated by a smaller 
momentum from Q2. On the contrary a thermal gradi- 
ent changes the number and the frequency of phonons 
coming from Qi and Q2 and we easily obtain: 



dt 



dijiLo) dv 



(3) 



°^ dx dx ^ 

showing that within this framework a combination of 
thermal and velocity gradients can transfer physical mo- 
mentum along the cantilever. 




FIG. 1: (Color online) a) Heated cantilever beam, of length L 
and thickness t. The temperature increases along the 2;— axis 
from Ta to Tb- Three thin sections Qk, at intervals of length 
are shown in gray, b) The cantilever oscillates in its first 
bending resonance mode. When it goes through a zero cross- 
ing, the local velocity is directed along the 2;— axis and its 
amplitude increases along the a;— axis, c) Detail of sections 
Qk- The density and the frequency of the phonons increase 
from left to right, proportionally to the local temperature. 
As phonons travel for a distance If, one half of the phonons 
starting from Q\ and Q3 reach the central section Q2- 



To evaluate the possible implications of this effect on 
the dynamics of a resonator, we have to put Eq. ^ on 
a more solid ground, in close analogy with the theory 
of transfer of momentum in the kinetic theory of gases 
[l6| . We consider the point P = {x,y,z) in a crystal, 
where the temperature is given by a scalar field T{x, y, z) 
and the velocity by a vector field v(x,?/, z). We assume 
that a steady state has been reached and evaluate the 
' ' balance of momentum in a small sphere A of radius dA 
centered at P. We place a polar coordinate system at 
(2) P, and consider a small volume da — sin 61 dr d6' d(/) at 



the point r = (r, 0, 0), say Q, in this coordinate sys- 
tem (see Fig. 1 of the supplementary material 18]). 
Both volumes are small in the sense that their charac- 
teristic size is much smaller than the MFP //. For each 
phonon branch s, in dtr are r7,s(TQ, k)d(T phonons with 
wave vector k, where Tg is the value of the tempera- 
ture at the point Q. Each phonon starts a free path cjlf 
times per second, where c is the group velocity, there- 
fore a number of — ^^^^j^^dtr dt phonons will leave the 
volume da in the time dt. Of these, the number in a 
cone of base 7r(dA)^ will have paths directed towards the 
small sphere if all directions are equally probable. That 
is, '^4^^2 "f^s (Tq , k) ^ dg dt phonons will leave headed for 
the sphere along r. Of these, a fraction only will suffer 
a collision within the sphere A: if we consider for in- 
stance only phonons directed exactly along — r, the frac- 
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tion e~^^~'^^^^^f will succeed in reaching the sphere at 
P, while the fraction e"'-'"^'^'^-'/'-'' will cross it without 
impact. As a result the number of collisions within A 
is |e~''/'-f2^ , where the factor | accounts for the re- 
duction of the length of the path on the borders of the 
sphere. Hence the number of phonons with wave vector 
k starting in da and colliding within the sphere A is: 

dN = l^(dA)3n,(TQ,k)^-^e-'-/'/dadt (4) 

The first term is the volume of the sphere A, that we 
approximate as a cubic volume element dV = dx dy dz to 
allow the integration over the crystal 17[. The collision 
rate at P is then: 



^=n.(TQ,k) 



-'■/'McrdF 



(5) 



The components of the total rate of the momentum en- 
tering the element dV are obtained by integrating the 
product of equation ([5]) and Here k varies over the 
full phonon spectrum and the volume da, placed at the 
source point Q, over the volume of the body Vb- 



dp 

dt 



= dV^ 



da 



Vb 
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(6) 



where the direction of propagation of the phonons has 
been approximated as ik — — ir- This equation can be 
evaluated after some preliminary consideration. First, 
the volume integral can be extended over all space, be- 
cause only a sphere of radius comparable with a few 
mean free path actually contributes to the volume in- 
tegral, due to the exponential term e~^^^f . Second, we 
note that ^iw^, n<;(r, k) = jj" c.ys(k)dr, where c^s(k) is 
the contribution of the mode (ks) to the specific heat 
C„ = Sks ''''''(^)- Then the sum over the phonon spec- 
trum is J2ks 7ks /J'ct,s(k)dr = ^{T)CydT, where 
7(^) = Sks 7ksC«s(k)/C^ is the overall Griineisen pa- 
rameter. In general this integral has a complex func- 
tional dependence on the temperature [l8| . as 7(T) may 
change in a complicated way. On the other hand at very 
low temperatures (below 10 K) we have /J"'' 7(T)C^dT ~ 
7o -^Tg, with By determined by the relation Cy ~ ByT^ 
and 7o = limr^o 7(7^)- At intermediate temperatures 
(30-70 K for silicon) this relation remains a good ap- 
proximation with 7o an effective Griineisen parameter 
|18| . We then expand vq and Tq to I order in r around 
the point P, and obtain for the components of ^: 



dpa By 



(7) 



VPa 



da 



pi-f 



WTp ■ Vvpc 



where a = x,y, z and all derivatives are evaluated at the 
point P. The details of the calculation, with a discussion 
of H order terms not shown here, can be found in the 
supplementary material [l^ . 

To complete the momentum budget we must consider 
the rate of momentum loss due to phonons starting a new 
free path in the element d^. For each direction {9, (/)) a 
number of dV ns(Tp, k)^ eded4> phonons will leave 
the element in the time dt, each carrying away the mo- 
mentum given by Eq. ([1} . It is straightforward to show 
that this contribution, when integrated over the phonon 
spectrum and all directions, compensates the first term 
(thermal equilibrium term) of Eq. ([7]). 

The second term (phonon pressure term) is propor- 
tional to the temperature gradient and describes the bal- 
ance of momentum along the direction of propagation 
ik. In fact, the distribution of phonons arriving from the 
hot part is shifted toward higher frequency, therefore the 
resultant momentum rate is directed against the temper- 
ature gradient. The force associated to this momentum 
rate is supported by the crystal through the elasticity 
modulus and induces a nonuniform expansion across the 
volume of the solid body. Hence our model extends to 
the case of a thermal gradient the well known connec- 
tion between phonon momentum and thermal expansion 

ail. 

The third term (damping term) is proportional to the 
gradients of the temperature and of the velocity fields, 
and can have striking effects in the case of an oscillat- 
ing crystal with a thermal gradient. In fact the rate of 
momentum transferred from the phonon gas to the crys- 
tal is equivalent to a force fp = If the point P in 
the time dt changes its position by dxp, the work done 
by this force is fp • dxp. This work can be easily eval- 
uated if the system is oscillating with angular frequency 
LJo . In this case the components of the displacement field 
are Ua{x, y, z, t) = Gufa{x, y, z) sin(a;ot), where G„ is the 
magnitude of the amplitude oscillation and fa(x, y, z) the 
displacement function. The velocity field and its gradient 
are respectively Va{x,y, z,t) — ujoGufa{x,y, z) cos{ujQt) 
and Vva = wqGuV/q cos(a;ot). If we consider that the 
force fp is much smaller than the elastic forces at the ori- 
gin of the oscillatory behavior, we can safely assume that 
it does not alter in a significant manner the displacement 
field Ua{x,y, z,t). Then the work done in an oscillation 
cycle is simply 6Lp = Jq^^'^° fp • vpdt, and this local 
contribution can be integrated over the volume of the 
crystal to obtain the total work: 
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^oGl J2 I (T' (VT ■ V/„) U) dV 



(8) 

This equation is the main result of this paper and must 
be carefully analyzed to understand the possible impli- 
cations on the dynamics of the resonator. 

We note first that the work Lph is proportional to the 
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squared amplitude of the oscillation, just like the energy 
losses due to internal friction in the resonator. On the 
other hand it can be positive or negative, depending on 
the interplay between velocity and temperature gradi- 
ents, therefore the superimposed heat flow can sustain 
or damp the oscillation of the body. To compare more 
easily this effect with the usual dissipative phenomena 
in resonating bodies, we define the loss angle due the 
phonon momentum transport as cjjph — "^"ly^i where 
W = \pG\u:l Sv /adV' is the total energy stored in 
the resonator. From Eq. ([5]) we have: 

where positive values of (j)ph mean a damping of the os- 
cillation. Besides this phononic contribution, the energy 
loss in a resonator results from the contribution of sev- 
eral dissipation mechanism, that collectively determine 
the background dissipation ISW^g and the background 
loss angle (f>bg — ^^''^ ■ The most common dissipa- 
tion mechanism are thermoelastic and structural loss in 
the resonator and energy loss through the support. In 
this respect, we note that silicon resonators based on 
cantilevered beams have demonstrated a background loss 
as low as 10~® at cryogenic temperatures, with a repro- 
ducibility of about 20% [ii]. 

For the cantilever of Fig. [1] the temperature profile is 
T{x) = (Tb —Ta)x/L + Ta and the displacement is fully 
described by the the component fz- We use the variable 

= x/L, ^ e (0, 1), and write the modal shapes as: 

MO = [sin(/3) - sinh(/3)] [sin(/30 - sinh(/30] + 
+ [cos(/3) + cosh(/3)] [cos(/30 - cosh(^0] (10) 

where j3 = 1.875 for the first resonant mode [13] ■ The 

frequency of the oscillation is obtained as = 

where Y and p are respectively Young modulus and den- 
sity of the crystal. In Fig. [2]we show the expected damp- 
ing for the lowest frequency mode of a silicon cantilever 
(y = 169 GPa, p = 2329 kg/m^) of thickness 100 ^m and 
lenght 50 mm. We use the average MFP from measure- 
ments of thermal conductivity accumulation distribution 
21 1, and 7(T) from thermoelastic measurements 22, 23| . 
In this case Eq. (0) cannot be used below 60 K, as the 
MFP would become larger than the thickness of the can- 
tilever and boundary scattering dominate the transport. 

Now we can trace out the main features of the theory 
and evaluate its experimental relevance. At first we note 
that, with a relative gradient of 10%, the magnitude of 
4>ph can be well above the typical background loss in this 
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systems [2J]. Second, according to the sign of 7(T) and 
of the gradient, the additional loss (pph can augment or 
reduce the overall damping of the resonator. Therefore 
with the application of a thermal gradient it is possible 
to reach a regime where the dynamics of the resonator 
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FIG. 2: (Color online) Expected damping due to the phonon 
momentum transport for the first bending resonance mode (55 
hz) of the silicon cantilever. The curve at 60 K is obtained in 
the low temperature approximation (Eq. ((9)1) with 70 = —0.8, 
while curves at 100 K and 150 K are evaluated by taking in 
account the temperature dependence of the Griineisen pa- 
rameter • A positive thermal gradient can be obtained by 
heating the free end of the beam. Then, depending on the 
base temperature, the additional loss (pph can augment or re- 
duce the background damping of the resonator. Typical back- 
ground loss for a silicon cantilever of this size is (f>hg — 10"''. 



is dominated by the phononic damping. In this case, 
the loss angle can be much higher than expected back- 
ground, spoiling off the performances of the device, or 
could be considerably reduced, possibly driving the sys- 
tem towards region of dynamic instability. From another 
perspective this effect could allow the design of micro- 
resonators exploiting the transport of phonon momen- 
tum as a power source, drawing energy from the thermal 
gradient to sustain their oscillatory motion. We point 
out that phonon momentum effects turn out to be ac- 
cessible in a macroscopic system thanks to the extreme 
energy sensitivity of low loss resonant systems, while a 
wider variety of mechanical effects can be observed at the 
nanoscale [5(. 

We have extended the phonon momentum concept to 
the case of an oscillating crystal, and within this frame- 
work we have evaluated the dynamic behavior of an ar- 
bitrarily shaped resonator with an applied thermal gra- 
dient. In case of a cantilever we have found a number of 
features that are at experimental reach and could have a 
significant role in the performances of future devices for 
precision measurements. For this reason we hope that 
these ideas provoke new experimental and theoretical ef- 
forts aimed toward understanding the effect of phonon 
momentum on dissipative processes in resonators. 
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This PDF file includes the following material: 

• Coordinate system (Figure 1) 

• Temperature integral: approximation for silicon (Figure 2) 

• Volume integral: first and second order approximation 

• Volume integral: medium temperature approximation for silicon 



COORDINATE SYSTEM 




FIG. 1: Coordinate system used in the volume integration. 
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TEMPERATURE INTEGRAL: APPROXIMATION FOR SILICON 



The sum over the phonon spectrum in Equation (6) of the main text reduces to: 



F{Tq) = V 7k. / " c,s(k) = / " i{T)C,dT 



(1) 



where ^{T) — X^fcs 7ksCi,s(k)/C„ is the overall Griineisen parameter and may change in a complicated way with the 
temperature. For instance in silicon ^{T) takes negative values between 19 K and 124 K. On the other hand the 
strong temperature dependence of the specific heat C^, can make negligible the low temperatures contribution to the 
integral. 

In the case of silicon at very low temperatures (below 10 K) we have: 



B 



Fvlt{Tq) =~ lo^Tl 



(2) 



with = 0.591 J/(m^K**) determined by the relation Cv ^ ByT^ and 70 ~ 0.4 [J. As shown in Figure [2j:, this 
relation remains a good approximation in the range (30-70 K) if we use an effective Griineisen parameter 70 ~ 0.8. 

In the medium temperature range (70-200) K, the temperature dependence of 7(T) changes the functional form of 
F(T), that now is well approximated by a quadratic functions: 



Fmt{Tq) — -Fo + FiTq + F2TQ 
with Fo = 2.85 x 10^, Fi = -6.17 x 10^ and F2 = 2.57 x 10^. 



(3) 
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FIG. 2: Evaluation of the temperature integral F{T) = j^"^ l{T) CvdT for the case of silicon, a) Temperature dependence of 
the overall Griineisen parameter 7(r) 0|; note that 7(T) is negative from 18 K to 124 K, where silicon has a negative coefficient 
of thermal expansion, b) Specific heat Cv of silicon Q]. c) Solid line: numerical evaluation of F{T). Hollow squares: points 
from the low temperature approximation Eq. ((5)) but with an effective Griineisen parameter 70 — —0.8. Hollow triangles: 
points from the medium temperature quadratic approximation Eq. ([3}. 
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VOLUME INTEGRAL: FIRST AND SECOND ORDER APPROXIMATION 



From Equation (6) of the main text, after the average over the phonon modes, we write the rate of momentum 
entering the volume element dV as: 

f -dv-jA_|^^d.F(r„)-l.-/^,(-.i.- + v«) (4) 

In the very low temperature case we have F{Tq) = Fvlt{Tq) —— 70-^ Tq, with 70 = limT^o7(2^)- As shown in 
the Figure [21 this relation remains valid at higher temperatures with 70 an effective Griineisen parameter. Therefore 
Equation ^ becomes: 



dp _ By 370 
dt ~ A AttPsC^' Jy^ 



dV f daT^le-'-/'^(-ca, + VQ) (5) 



With reference to the coordinate system shown in Figure [1] we expand vg and Tq up to II order in r: 

1 2 

VQa = Vpa + r • Vvpa + - (r • V) Vpa (6) 

^ Tj, + r • Vr| + ^ (i" • C^) 

where the derivatives are evaluated at the point P, corresponding to r = 0. The volume element is da = sin OdrdOdcj), 
therefore the volume integral for the component a, with a = x, y, z, is: 

r-TT p2n / 1 \ / 1 



poo pTT pZTT / 1 

/ dr de / d(f, smde~''^^f [T^ + r- \7T^ + 7: {r ■ V)^ 
Jo Jo Jo \ 2 



-Cir ■ ia + Vpa + r ■ Vvpa + 2 ('^ ' '"Pa ) (8) 



where we have extended the volume integral over all space, because only a sphere of radius comparable with a 
few mean free path actually contributes to the volume integral, due to the exponential term e~'^l^i . Here = 
sin Q cos ia; + sin B sin (/) iy + cos d and r = ri,. 

We note that, as derivatives are evaluated at a fixed point, each term is the projection of a constant vector along 
ir- Then, in the angular average: 



• a single projection averages to zero due to the angular oscillating term: 

p27T pTT 



/ d(f> sin^de* [ir • ia] 
Jo Jo 

/ d0 / smed0[ir-Wh{x,y,z)]^O 
Jo Jo 



(9) 
(10) 



where fi{x, y, z) is a scalar field and the gradient is evaluated at r = 0. 

• The product of two projection gives a finite value: 

J^^ d<j> sin e de [ir ■ io] [i. • V/i {x, y,z)] = ^n ^^^^g^^'^^ 

/ d(j) sm9d9[ir ■'^fi{x,y,z)][ir ■Vf2{x,y,z)] = -TTWfi{x,y,z) ■Vf2{x,y,z) 
Jo Jo -J 

/ dcf, f sin 9 d9 {ir • V) Vi y.z) = \TT VVi V, z) 
Jq Jo •J 

where f2{x, y, z) is another scalar field and again the gradients is evaluated at r = 0. 

• The product of three projection gives again a null value: 



(11) 
(12) 

(13) 



/ del) sm9d9[ir-ia,] [v-Vf h{x,y,z) =0 (14) 
Jo Jo '- -' 
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As an example, we evaluate Equation (ITT|) in the case a = x. Now [i^ ■ !«] = s'm9 cos cj), therefore 

/ d(f) [ sin6'd6' I sm^cos^^^ + sin6'sm(/<^^ + cos^^^ ) sm0cos0 = 
Jo Jo \ ox dy dz ) 

d(j)J^ sin^de'^sin^eicos^,/)^ 



9h .... r.„3..._4 9A 



, cos d(/) / sin"^ d0 = - tt ■ 
The integration over the radius r is much simpler: 



e 





(15a) 




(15b) 


""'^f dr = 21) 


(15c) 



The rate of momentum entering the volume element dV^ is then: 

%i ^ dy 70^ - 4T| ^ + ^n^-l VTp ■ Vvp^ + 12^ vp^i^Tpf + ^TJ. V2^;p„') (16) 

at 4 \ IfC oa c c c J 

In respect to the I order expansion in r leading to the Equation (7) of the main text, now all contributions up to the 
first order in the mean free path If are considered. Further terms of the expansions ^ and ([7]) would add negligible 
contributions proportional to I'j. We have already discussed the first three elements of this equation, respectively 
the "thermal equilibrium", "phonon pressure" and "damping" term, and address now the relevance of the additional 
terms. 



• The fourth term is a correction to the damping term, as it depends both on the temperature gradient and on 
the local velocity of the crystal. In the case of a silicon cantilever discussed in the main text, the correction 
is less than 10% over the considered range of relative gradient. We note that to derive this term we used the 
relation: 

(T^) = 3T}{VTp f + TlV'^Tp ^ 3 T^{\7Tp f (17) 
where S/^Tp = because Tp is a solution of the stationary heat equation with no source term. 

• The fifth term is a correction to the thermal equilibrium term, therefore the associated loss angle concurs to 
determine the background loss. It can augment or reduce the overall damping of the resonator, depending 
on the modal shape and on the Griineisen parameter 7 of the lattice. Lacking a the dependence from the 
gradient, we believe that it is not reasonable to single out this phononic contribution among the other dissipative 
mechanism. As a consistency check, we only note that for the silicon cantilever discussed in the main text the 
expected contribution is a few units xlO~^, not far from the values measured in similar systems [4|. This term 
is consistent also with the extremely low loss factor measured in resonators exploiting torsional modes, such as 
the Double Paddle Oscillator [5]. In fact in this case the expected damping from phonon transport is orders of 
magnitude smaller due to the symmetry of the modal shapes. 
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VOLUME INTEGRAL: MEDIUM TEMPERATURE APPROXIMATION FOR SILICON 



From Equation (j4|, in the medium temperature case, we write the rate of momentum entering the volume element 
dV as: 

^ = 4^ 1 d.^...(T,)-L e-^/'. (-ci. + V,) (18) 

with Fmt{Tq) given by Equation In this case it is straightforward to show that, to the first order in Zj, the 
contributions to the rate of momentum entering the volume element dV are: 

• Thermal equilibrium term 

dV {^-^ + ^V\p„) (Fo + FiTp + F2TI) (19) 

• Phonon pressure term 

- dF ^ (Fi + 2F2TP) (20) 

• Damping term 

21 t 21 f 

AV (Fi + 2F2TP) VTp ■ V«p„ + AV -^F2Vp^{VTpf (21) 
c c 

These local contributions were integrated over the time and over the volume of the cantilevered beam to obtain the 
total loss angle at 100 K and 150 K reported in Figure 2 of the main text. 
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